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Abstract. New calculations of the Galactic contamina- 
tion over microwave background radiation anisotropies are 
carried out. On one hand, when a frequency-dependent 
contrast of molecular clouds with respect to the Galactic 
background of the diffuse interstellar medium is taken into 
account, the anisotropic amplitude produced by Galactic 
dust is increased with respect to previous calculations and 
this is of the same order as that of the data from the ob- 
servations. On the other hand, if we take into account 
rotational dust emission, for instance, a frequency inde- 
pendence of anisotropies in the microwave range may be 
obtained. 

This leads to the possibility that under some partic- 
ular, but not impossible, conditions all the microwave 
background radiation anisotropies may be due to Galac- 
tic foregrounds rather than cosmological in origin. More- 
over, a suspected coincidence between the typical angular 
sizes of the microwave background radiation anisotropies 
and those of nearby molecular clouds makes more plau- 
sible the hypothesis of a purely Galactic origin for these 
anisotropies. It is also argued that the correlation among 
structures at different frequencies, the comparison of the 
power spectrum at different frequencies and the galactic 
latitude dependence of the anisotropies are not yet proofs 
in favour of either a cosmological or Galactic origin. 
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1. Introduction 

Various observations of the cosmic microwave back- 
ground radiation (CMBR) have led to the claim that 
it is anisotropic (see reviews Readhead & Lawrence 
1992; White, Scott & Silk 1994). Measurement of these 
anisotropies provides information on structural formation, 
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inflation, quantum gravity, topological defects (strings, 
etc.), dark matter type and abundance, the determination 
of Hq, il, A, the geometry and dynamics of the Universe, 
the thermal history of the Universe at the recombination 
epoch, etc. A knowledge of these anisotropies is very im- 
portant for discriminating among different cosmological 
models as well as for measuring certain related parame- 
ters. 

Unfortunately, this information is contaminated by 
several effects which are unrelated to cosmology, and 
which are either extragalactic or Galactic in origin. Ex- 
tragalactic sources of such contamination may be gravita- 
tional scattering due to inhomogeneities in the matter dis- 
tribution of superclusters (Fukushige, Makino & Ebisuzaki 
1994; Suginohara, Suginohara& Spergel 1998), inhomoge- 
neous reionization (Knox, Scoccimarro & Dodelson 1998), 
possible radiative decay of massive tau neutrinos (Cele- 
bonovic, Samurovic & Circovic 1997), the inverse Comp- 
ton effect, temporal variations of the gravitational poten- 
tial, etc. This paper will not analyse such extragalactic 
sources of anisotropy but will concentrate instead on those 
coming from our own Galaxy, particularly at large angular 
separations (9 > 2°). 

A summary of possible contributions from our own 
Galaxy to the anisotropies has been studied, for instance, 
by Bennett et al. (1992) and has led to the conclusion 
that there are three possible emission mechanisms: syn- 
chrotron, free-free and dust emission. The calculation of 
anisotropies originating in our Galaxy was carried out 
on the basis of observations from the COBE-DIRBE, 
which were extrapolated to other frequencies, and it was 
concluded that Galactic contamination is negligible. The 
paper of Bennett et al. was optimistic. However, not 
all authors think the same. Masi et al. (1990) pointed 
out that dust may provide an important contribution 
to the anisotropies in the microwave region. Banday & 
Wolfendale (1991) predicted that dust contamination is 
quite important except in some regions away from the cir- 
rus. Hence, whether or not the dust contribution is negli- 
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gible seemed to be a point of contention by the beginning 
of nineties. 

Recently, clear evidence of correlations between far- 
infrared maps (which trace Galactic dust) and microwaves 
has been presented (Kogut et al. 1996a, de Oliveira-Costa 
et al. 1997; de Oliveira-Costa et al. 1998). In the light 
of these works it is difficult to deny the influence of the 
Galaxy on the microwave background. Even though the 
source of the contamination remains unknown, the corre- 
lation between the Galaxy and COBE-DMR data remains 
an observational fact. Another recent paper written by 
Pando, Valls-Gabaud & Fang (1998) points out the non- 
Gaussian nature of the CMBR over large angular scales 
of anisotropics against standard cosmological predicitions, 
and this might be due to an excess of foreground contam- 
ination. 

In this paper, an explanation is offered for these last 
observational facts. Two new elements are added to pre- 
vious works for the calculation of Galactic contamination: 
a frequency-dependent contrast of molecular clouds with 
respect to the Galactic background of the diffuse interstel- 
lar medium; and the existence of a recently noticed source 
of emission: rotational emission by dust grains (Draine & 
Lazarian 1998a). These lead to a non- negligible and pos- 
sibly unique source of anisotropies due to our Galaxy. 

2. How large is the Galactic contamination? 

The existence of flux anisotropies due to Galactic clouds is 
undeniable. In the far infrared, the presence of the molec- 
ular clouds is associated with "infrared cirrus". There is 
a good correlation between infrared cirrus and atomic hy- 
drogen gas (Boulanger, Baud & van Albada 1985; Bur- 
ton & Deul 1987; Boulanger & Perault 1988; Schlegel, 
Finkbeiner & Davis 1998) and the fluctuations of the hy- 
drogen column densities within the nearest 75 to 100 pc 
of the Sun are occasionally about 5 x 10 19 cm~ 2 (Frisch & 
York 1986), so this implies that fluctuations due to nearby 
clouds must be detected to some extent in the microwave 
region. 

Most of the clouds are confined in the Galactic plane 
but there are few which can be observed at high galactic 
latitudes, and which are close to us (Blitz 1991). Those re- 
sponsible for large angle correlations would need to have a 
large size (i.e. giant molecular clouds). So, some microwave 
background radiation anisotropies (MBRAs) may be due 
to inhomogeneities in the density distribution of the local 
interstellar medium. 

Dust associated with these giant clouds exists in three 
different stages associated with molecular, neutral and 
ionized hydrogen (Sodroski et al. 1997), whose total emis- 
sion gives a continuum spectrum proportional to the col- 
umn density in the line of sight. Column density fluc- 
tuations lead to intensity fluctuations. Dust anisotropies 
would be due to its non-homogeneous distribution in 



molecular clouds or its equivalent for other types of Galac- 
tic emission. 

Dark matter around spiral galaxies in the form of 
cold gas, essentially in molecular form and rotationally 
supported, was purposed by Pfenninger, Combes & Mar- 
tinet (1994) and Pfeninger & Combes (1994); this cold 
gas in molecular form was also held to contribute at mi- 
crowave wavelengths (Schaefer 1994; Schaefer 1996). The 
gas would be at temperatures close to 3 K and composed of 
small clumpuscules of size 30 AU. It is claimed (Combes & 
Pfenniger 1997; Schaefer 1996) that part of the radiation 
observed by COBE-FIKAS is cold gas instead of dust. It 
might be in the outer disc and difficult, although not im- 
possible, to detect. Indeed, Combes & Pfenninger (1997) 
propose techniques for detecting it. Though this cold gas 
may be a source of anisotropies, it would be only for very 
small scales since the emission is high enough just for very 
high densities (n(7?2) ~ 1 18 cm -3 ) inside the small clum- 
puscules. Thence, cold molecular gas emission will not be 
further considered. 

In following subsections, the expected emission from 
our Galaxy will be compared with observations in the mi- 
crowave region. 

2.1. MBRAs amplitude 

The two-point correlation function is: 

(T*T*)(6) = [ dtt [ dQ'T*(Q)T*(fl') (1) 

J J\o'-a\=e 

, where T* is the .antenna temperatureS, once the average 
flux is subtracted!], the results are those shown in Fig. |^ 
and m 

In order to calculate the effects of the Galaxy on 
MBRAs, extrapolation of the measured two-point corre- 
lation functions (Fig. |^) at 1.42 GHz due to synchrotron 
emission and at 240/im due to continuum dust emission 
must be carried out. Free- free emission is not considered 
here in order to simplify the calculations and because we 
have no information concerning it (Smoot 1998). In any 

1 Antenna temperatures, T, were obtained from the intensity 
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maps at 1.42 GHz and 240/im multiplied by the factor If a , 
where c is the light speed and k is the Boltzmann constant. 

2 A background depending on Galactic coordinates is re- 
moved to eliminate the flux variation due to its smooth gra- 
dient. Therefore, we measure only the correlations due to flux 
anisotropies with regard this average background. I have used 
bidimensional "spline3" functions of order 3 in both coordi- 
nates and with crossed terms to fit the background in the 
selected regions regions. This was done by means of the IM- 
SURFIT task of IRAF. These regions are selected in off-plane 
regions (|6| > 20°). See the map at 240/im with the background 
subtracted in Figure |l]. For 1.42 GHz, since the whole sky was 
not covered, only available regions with \b\ > 20° were used 
(S > —19°). The zodiacal component is not removed at 240 
fim as it is negligible at these frequencies (Reach et al. 1995). 
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Fig. 1. The COBE-DIRBE flux map at 240 /im after subtrac- 
tion of a smooth component of the Galactic background as 
explained in the text. The contours are from -20 to +20 MJy 
sr _1 in steps of 4 MJy sr _1 . 
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Fig. 2. Angular correlation function for different surveys 
at \b\ > 20°: a) 1.42 GHz (the Max Planck Institute for 
Radio Astronomy, Reich 1982; Reich & Reich 1986), the 
dashed line is a fit to 0.0663 - O.OO3820(deg.) K 2 ; b) 240 fim 
(COR&DIRBE, Boggess et al. 1992), the dashed line is a fit 
to 2.13 x HT^deg.)- 1 - 20 K 2 . 

case, if it is proven than dust or synchrotron emission or 
both are high enough then adding free- free emission would 
increase Galactic anisotropies, never decrease. 

At any frequency in the microwave range, two-point 
correlation functions must be evaluated by extrapolating 
both effects independently and summing them up accord- 
ing to: 

(T{v)T{v)) = ((T synchr (z/) + T dust (v)) 

X(7synchr(^) + 7d us t ("))) = {T sync Yn{ v )Tsynchr (^)) 
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Fig. 3. Cross correlation function for two surveys: 1.42 GHz 
(Max Planck Institute for Radio Astronomy, Reich 1982; Reich 
& Reich 1986) and 240 /im (COB.&DIRBE, Boggess et al. 
1992). 



+ (Tdust(^)Tdust(^)} + 2(T synch r(^)T , dust (f)). (2) 

Each of these terms may be calculated by means of an 
extrapolation of previously measured correlations. I will 
calculate the contribution due to the local Galactic flux 
fluctuations with respect the average flux, i.e. < T*T* >, 
and not that due to the smooth variation of the depen- 
dence of this average flux with Galactic coordinates. Usu- 
ally, the extrapolation is carried out by multiplying each 
of the pixels of the map at Vq by the mean amplitude 
variation of the flux {(T{y) /T(vq))) and calculating the 
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WRONG T 
EXTRAPOLATION 90 ' 

(the one used by the rest of the authors) ^ 





EXTRAPOLATION 
OF THIS PAPER 

(it takes into account that the fluctuations 
and the average emission are not equally 
proportional for different frequencies) 



0.1 



Fig. 4. Example of extrapolations for Galactic dust emission. 
The main difference between the extrapolation of expression 
(^) and those derived in work prior to this paper is that the 
latter do not take into account the frequency dependence of 
the proportionality between the amplitude of the fluctuation 
and that of the average flux. 

correlations of the new map at frequency v. This is equiv- 
alent to multiply the two-point correlation function by a 
factor {T{i/)/T{vq)) 2 . I will go one step further and will 
take into account that the local fluctuations do not vary 
proportionally to the average flux (Fig. ||); therefore, the 
mean amplitude of the two-point correlation function is 
affected by a second factor. As we will see later, this vari- 
ation is due to different effective dust temperatures be- 
tween the Galactic clouds which produce the fluctuations 
and the Galactic medium that contributes to the average 
flux. The shape of two-point correlation f unct ion is also 
expected to vary, as will be explained in §3^3. However, 
since we want to estimate only the order of magnitude of 
the correlation amplitude in the microwave region, we will 
not take this effect into account at this stage. The effects 
of the shape variation of the t wo-point correlation func- 
tion will be analysed in § 3.3.1 . Hence, extrapolations are 
calculated according to 



(Tgynchr ( V) T sync hr \v) ) (6) 



/ AT synchr M \ 
V T 3ynchr M J 

' AT Bynchr (1.42 GHz) 
T aynchr (1.42 GHz) 



-^synchr(^0 



T synchr (1.42 GHz) 



; <r; ynchr (i.42 GHz)r; ynchI .(i.42 gh z )>(0), 



(3) 



Jdust(^) \ 2 

r dust (240 H/ 



V T duat H 



/ AT dust (240 M m) 
\ T dust (240 fj,m) 



x(T d * ust (240 M m)T| ust (240 M m))(0), 

(Tsynchr (^)Tdust (v)) {&) 

— / ( ^synchr(^) ?dust(^) 

~ \ Vr Bynchr (1.42 GHz) T dust (240 /xm) 



(4) 



/ AT 3ynchr H 



( AT dust Q) \ 
V T duBt («/) ) 



( AT 3ynchr (1.42 GHz' 
, V T synchr (1.42 GHz) 



( AT du3t (240 M m 
V ^du St (240 /mi) 




x(T s ; nchr (1.42 GHz)T d * ust (240 fim))(9), 



(5) 



so we have to multiply three factors for any of the corre- 
lations: the first factors ((T(v)/T{vq))) take into account 
the variation of the mean Galactic background, the sec- 
ond factors ({^-{v) / ^{vq))) for the variation of relative 
fluctuations at different frequencies and the third factors 
(< T*T* > (vq,8)) are respectively the correlations mea- 
sured in Fig. a), b) and |[ The three factors can be 
separated because of the assumption of the independence 
of 9 for the first and the second factor. 



First factors: The antenna temperature due to 
chrotron emission is 



syn- 



Tsynchr^, |6| > 20°) = 0.6 



K, 



(6) 



. 1420 MHz 

where (3 « —3 (Davies, Watson & Gutierrez 1996). See 
the contribution in Figure [|. 

The dust temperature has two contributions: thermal 
emission, which is detected at far infrared wavelengths 
(this is the emission detected at 240 /jm); and the 
microwave rotational emission due to small spinning 
grains (Draine & Lazarian 1998a). The first contribu- 
tion is widely accepted while the second is still usually 
not included in the calculation of the Galactic contri- 
bution to the anisotropics in the microwave range. 
The mean thermal dust emission is given (Fig. ||) by 
Reach et al. (1995): 

7d us t>, l&l > 20°) = ( - Y [1.74 x 10 -5 

dusty ,n ) ^ 900 GH J L 

xB(v,T = 17.72 K) + 1.23 x 10~ 4 B(^, T = 6.75 K)] 



2kv 2 



K. 



(7) 



where B is the blackbody radiation intensity, c the 

2 

light speed, k the Boltzmann constant and ^f^y is the 
factor for converting intensities into antenna tempera- 
tures. This is an average law for \b\ > 20°; as shown by 
Reach et al. (1995), there may be important deviations 
from this law for some regions. 
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The second term of eq. (Jfy may be due to very cold 
clouds; their average temperature would be 6.75 K 
if dust emission were proportional to v a B{u, T) with 
a = 2, although this would be somewhat different if 
a had another value. Expression (fy is just a fit to 
the observational COBE-FIRAS data, not a model, 
and their parameters may not have any direct phys- 
ical meaning. It was argued by Lagache et al. (1998) 
that an alternative hypothesis to the existence of these 
very cold clouds is a Cosmic Far Infrared Background 
due to distant galaxies of 26 /iK in the range between 
400 /zm and 1000 ^m (Puget et al. 1996)1. Regardless 
whether or not this hypothesis is true, for our purposes 
the extragalactic emission can be added to the aver- 
age flux of the Galaxy. The important thing, as will be 
seen later, is that the fluctuations of the flux are due 
to local clouds, and that the average temperature of 
the sources which are origin of the fluctuations is less 
than that corresponding to the diffuse emission of the 
Galaxy (plus any extragalactic background). 
The rotational emission is the predicted model by 
Draine & Lazarian (1998a; preferred model: A), Fig. 
U in this paper, taking into account a hydrogen col- 
umn density N(H) = 3.9 x 10 20 [cosec(6) -0.17], which 
includes neutral and ionized hydrogen (Heiles 1976; 
Reynolds 1991), and averaged over |6| > 20°. 
Second factors: I stress that, as far as I know, all au- 
thors who make the extrapolation of correlations set 
the second factor equal to unity for all frequencies, 
i.e. the relative fluctuations for dust or synchrotron do 
not depend on the frequency (for example, Banday & 
Wolfendale 1991; Bennett et al. 1992; Guarini, Mel- 
chiorri & Melchiorri 1995; Reach et al. 1995; Femema 
et al. 1998). This is a very poor approximation for 



3 This extragalactic background is obtained after subtracting 
the CMBR, zodiacal emission and Galactic dust emission from 
COBE-FIRAS data, making an extrapolation of the HI column 
density data with an emission law v 2 B(u,T — 17.5K). Obvi- 
ously, if it is assumed at first that there are not very cold clouds 
then the conclusion must be that there are no very cold clouds, 
and this is what these authors find, attributing the remnant 
emission to an extragalactic source. This remnant emission is 
quite isotropic and led Puget et al. (1996) and Boulanger et 
al. (1996) to think that this component does not belong to the 
Galactic disc; however, it is difficult to tell what really happens 
because of low signal-to-noise at high Galactic latitudes and 
neither does Galactic cirrus have a clear position dependence 
at high latitudes (see §3.4); thus, I think these local sources 
are not excluded from containing very cold clouds. Moreover, 
the extrapolation of HI column density to far infrared emis- 
sion is also dangerous since the correlation is not rigorous (see 
3). HI is a dust tracer but not a perfect one and it may be 



not appropriate at all to find out information about the origin 
of some small remnants of COB-&FIRAS which do not corre- 
late with it. From all these considerations, I think it is risky 
to place the origin of some very far infrared radiation (400 to 
1000 micron) outside our Galaxy. 
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Fig. 5. Mean Galactic background emissions for |6| > 20° 



dust, especially when the extrapolation is to one or- 
der of magnitude for lower frequencies, as it is in this 
case. I think this poor approximation is the main rea- 
son why dust anisotropies were considered negligible 
in the past. 

The reason why I think this is a bad approximation for 
dust is that the mean dust emission intensity follows 
a different continuum spectrum from that of the local 
clouds (infrared cirrus) providing fluctuations at scales 
of several degrees. When the dust emission is modelled 
by I oc v a B{v, T) (with B the blackbody radiation), a 
varies from place to place (Tegmark 1998): it is smaller 
for molecular clouds (Mathis 1990; Schloerb, Snell & 
Schwartz 1987; Meinhold et al. 1993) or higher lati- 
tudes (Banday & Wolfendale 1991; Reach et al. 1995) 
than that for the mean intensity (Matsumoto et al. 
1988; Wright et al. 1991; Reach et al. 1995). Moreover, 
the mean intensity includes emission from star-forming 
regions and mass-loss stars in which the dust is hotter 
than the diffuse or clouds dust (Draine 1994). The a 
spectral index depends on the grain properties (Wright 
1993) varying from 1 to 3.5 (Banday & Wolfendale 
1991). It is Rs 2 for graphite, w 1.5 for silicates and « 1 
for layer-lattice materials, such as amorphous carbon 
(Draine & Lee 1984). The difference is also related to 
the different distribution of dust-grain sizes in the dif- 
fuse interstellar medium and in molecular clouds. The 
rate of small grains is larger in the diffuse interstellar 
medium than in the clouds and the temperatures are 
different implying that the peak emission from clouds 
occurs at different wavelength as compared with that 
of the diffuse medium (Greenberg & Li 1996). 
These second factors are unknown because the 
anisotropies are not sufficiently explored for wave- 
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lengths longer than 240 fxxsx. In order to give a rough 
estimate, a simplifying assumption can be made: the 
Galactic contribution is the sum of the local clouds — 
infrared cirrus — (T c ) plus a diffuse Galactic back- 
ground (Id), where all the anisotropies at angular 
scales of several degrees are due to anisotropies in the 
local clouds distribution (AT d = 0). Hence, 



Tdlist — T, 



dust.d 



T. 



dust,c 



( AT du3t , c (i/) \ 
V T du3t ,<» ) 



1 



(8) 



(9) 



The numerator of this expression is constant for all 
frequencies because the fluctuations of the tempera- 
ture due to clouds are proportional to their density 
fluctuations for any frequency. However, the denomi- 
nator is not constant. The variations of the rate Td/T c 
as a function of the frequency must be estimated to 
evaluate the second factor. 

Dust thermal emission. For a wavelength of 240 
jitm the rate can be directly measured from the 
COBE-DIBE maps at |6| > 20°, taking T as the 
total radiation, which includes diffuse and clouds 
emission; and T* as the radiation once the back- 
ground is subtracted (as explained above), which 
is T* = T c - (T c ). If we assume (|AT C |) - (T c ), 
i.e. the typical fluctuations in the cloud flux is as 
large as the average of the flux, which is justified 
in a distribution of few single clouds over \b\ > 20° 
with many regions having a negligible contribution, 
then 

r dust ,d(240 M m) J|b|>20° ^r(Q) 
+ T dusM (240 A1 m) ~ Jj 6|>20Q <ffl|T*(fi)| ) 
The outcome implies that 30% of the radiation at 
240 ^m at |b| > 20° comes from clouds and 70% 
from the diffuse interstellar medium. This is just 
an estimation, but it is not a very bad approach. A 
comparison could be made with the rate from the 
model by Cox, Krugel & Mezger (1986) which gives 
a value even greater: ~ 4.5 in the Galactic disc. 
Clouds are colder than the diffuse medium: between 
6 and 15 K for molecular clouds and around 16 
K for diffuse medium according to Greenberg & 
Li (1996). Not all the infrared excess can be due 
to molecular clouds, but all of these are generally 
colder than the diffuse medium of dust associated 
with gas. Indeed, infrared-excess clouds are peaks 
of column density of dust probably associated with 
molecular gas of colder temperatures than the rest 
of the dust (Reach, Wall & Odegard 1998, Lagache 
et al. 1998). Hence, the increase in emission from 
clouds is 10 or 15 times greater than from the inter- 
stellar medium when we compare 240/mi and the 
microwave region, v < 100 GHz, and it is domi- 
nated by clouds emission in the microwave region 



(see, for example, Fig. 14 of Beichman 1987; or Fig. 
1 of Cox, Krugel & Mezger 1986) and 
1 T dustid (v < 100 GHz) _ 

T dust , c {v < 100 GHz) ' 1 ' 

Thus, the second factor for thermal radiation for 
these microwave frequencies is 

/ AT du3t (i/<100 GHz) \ \ 

[TC^wj (12) 

/ AT du3t (240 M m) \ I V ; 

^ T du3t (240 Mm) J / 

Dust rotational emission. A factor must also be 
taken into account for the rotational emission of 
the dust. First, an error in the function given by 
Draine & Lazarian (1998a) may be included in this 
factor since the paper uses different parameters 
that are poorly known. Secondly, the rate between 
small particles and large particles is lower for clouds 
than for the diffuse interstellar medium (Greenberg 
& Li 1996) because there are processes that de- 
stroy small grains (Puget, Leger & Boulanger 1985; 
Puget & Leger 1989), and this would decrease the 
emission. Thirdly, the density in the clouds is also 
higher, thereby providing higher emission. There is 
not enough data in the literature to carry out ac- 
curate calculations. In any case, my intention is to 
show the order of magnitude of the Galactic con- 
tribution rather than make accurate predictions, 
so I think it is enough to take the dust rotational 
emission shape of Fig. ||, divided by Tdust (240/zm), 
and multiply this by some free parameter /. The 
range of this parameter will be between and 1, 
so the intensity of the rotational emission is not 
overestimated. It will be conservatively estimated 
in the range of possible values compatible with our 
knowledge of clouds and rotational emission. The 
factor has to be less than unity because the values 
of the anisotropies predicted with a factor 1 exceed 
the total observational values for anisotropies mea- 
sured by the COBE-DMR and other surveys. 
Summing up, the second factor in the microwave region 
is taken to be 1 for synchrotron emission, 3.3 for ther- 
mal dust emission and an unknown factor < / < 1 
for rotational emission. The calculation of the second 
factor corresponding to the total dust emission is cal- 
culated for each frequency as a weighted average of the 
thermal and rotational emissions whose weights are the 
respective intensities of the first factors. 

So, with all this information the question to solve is 
how large are anisotropies due to the Galaxy at microwave 
frequencies. The answer is in the eq. (^|). For example, at 
90 GHz, the results are those shown in Fig. ^. Thermal 
emission of dust is the predominant one in this range. For 
this frequency, rotational emission with / < 1 does not 
contribute significantly, so the exact value of the factor 
/ may be ignored. Synchrotron effects are also negligible 
since the relative fluctuations in Fig. Ha) are quite low. 
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90 GHz 

|b|>20° 




degrees 



Fig. 6. Predicted Galactic two-point correlations at 90 GHz, 
\b\ > 20°, according to eq. (|) compared with the COBE-DMR 
data (Hinshaw et al. 1996). 

As can be observed in Fig. EL the amplitude of the 
anisotropy is of the same order as that observed by COBE- 
DMR. This casts doubts on the origin of MBRAs as . as 
will be discussed in section || they may be totally of Galac- 
tic origin. 

The results from ULISSE experiment (Bernardis et 
al. 1992) or other previous experiments (Melchiorri et al. 
1981) at 6° scale show that CMBR anisotropies at 90-870 
GHz or 410 GHz respectively are less than 35 or 40 /iK 
at some regions at high Galactic latitudes, after Galactic 
emission subtraction. However, these authors use a Galac- 
tic emission model which does not take into account the 
cold molecular clouds discussed in this paper. From eq. 
(|), averaged over \b\ > 20°, V /(TT)(6°)(870 GHz) ~ 
50/iAT (taking the "second factor" as / ~ 1 since 
cold/diffuse clouds flux ratio varies negligibily at these 
frequencies) and y/{TT) (6°) (90 GHz) - \9fiK (/ ~ 3.3). 
Fluctuations for the intermediate frequencies cannot be 
calculated since the value of / (second factor) as a func- 
tion of the frequency is unknown. Nevertheless, it can 
be seen that the order of magnitude of expected CMBR 
anisotropies and those from the Galaxy as calculated in 
this paper are nearly the same. 

3. Is it possible that MBRAs are purely Galactic 
in origin rather than cosmological? 

Maybe under some particular, but not impossible, condi- 
tions all the microwave background radiation anisotropies 
be due to Galactic foregrounds rather than cosmological 
in origin. Some arguments could be given against this and 
I will discuss them in next subsections. 



3.1. Frequency dependence of the MBRAs amplitude 

The main argument in favour of a cosmological origin 
of the MBRAs is that these do not depend on fre- 
quency. Strictly speaking, observations point out that the 
anisotropies at 53 and 90 GHz are nearly the same, but 
that they are higher at 31 GHz (by a factor 2 or so in dif- 
ference, Hinshaw et al. 1996). This favours a cosmological 
origin of the MBRAs but not to the total exclusion of a 
Galactic origin. 

The main argument against Galactic anisotropies with 
no frequency dependence in the range between 50 and 
90 GHz has been that any emission — dust, free- free or 
synchrotron — -gives a dependence proportional to v a with 
a 0. However, we should consider that there may be 
some range of frequencies of transition between a > — 
thermal dust emission — and a < — other kinds of emis- 
sion. The introduction of a quite potent rotational emis- 
sion (see Fig. ^J) is not normally taken into account and 
it is another important element, which has a negative a 
for v > 20 GHz. In the intermediate range between both 
regimes, the amplitude of the anisotropies follows approx- 
imately a constant dependence with frequency, i.e. a sa 0. 

Is it possible to get from Galactic sources of 
anisotropies the same two-point correlation function for 
53 and 90 GHz as observed by COBE1 The answer is 
yes: when / w 0.8 (see Fig. |?]). As observed in the fig- 
ure, there is a coincidence in the multiplication of the 
first and second factors contained in eq. (^) for 53 and 
90 GHz, so the amplitude of the correlation will be the 
same. Synchrotron contribution is low. There is also a 
dip around 60 GHz which separates slightly from a flat 
spectrum but this range has not yet been observed ac- 
curately. With the present model of rotational emission 
(model A of Draine & Lazarian 1998a), there is an excess 
of anisotropies for 31.5 GHz with respect the COBE-DMK 
observations. This, however, might be reduced when we 
choose other parameters for the dust rotational emission 
predictions by Draine & Lazarian (1998a). The aim here 
has been to show that it is perfectly possible to explain 
the observational anisotropies at microwave frequencies in 
terms of Galactic clouds rather than constructing an accu- 
rate model. The rotational emission is not well known and 
there is a lot of work that remains to be done to get exact 
results. Up to now, only rough estimates can be carried 
out and these indicate that, contrary to what was believed 
in the past, there exists the possibility that the totality of 
the MBRAs are Galactic. 

The main objection to this argument might be that 
it would be an enormous coincidence that / should be 
0.8 or, equivalently, in any rotational-emission prediction, 
that the range of nearly constant a should be between 50 
and 90 GHz. 

The predicted existence of rotational emission much 
more intense than free-free emission (Draine & Lazar- 
ian 1998a) is controversial and not yet proven. Its exis- 
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One remarkable feature of MBRAs that rouses suspicion 
about their relationship to our Galaxy is the coincidence of 
the typical angular size of their structures with the typical 
angular size of nearby clouds. These structures have an 
appearance very similar to the clouds observed in other 
frequencies. 

As an example, compare Fig. 4 of Gutierrez et al. 
(1997, reproduced in Fig. ^ c)), showing structures ob- 
served by the Tenerife Experiment, and Fig. || a), b), d), 
e) for other frequencies. The second differences are evalu- 
ated according to 

Second diff.(a,(S) = T(a,6) - -[T(a + 8.1°/ cos(5), 5) 



+T(a - 8.T7 cos(<5), S)\, 13 



Fig. 7. Multiplication of first and second factors of eq. (|2j) 
for dust anisotropics when the parameter / is equal to 0.8. 
Normalized such that AT(240/wn) = T(240^m). 

fence is used in this paper as a plausible candidate for a 
non-thermal contribution of emission correlated with dust. 
However, it is not essential to justify anisotropies at 53 
GHz. Other mechanisms could be used instead. Free-free 
is an option, though there are problems in justifying the 
necessary flux for the total anisotropies. Another mech- 
anism is the thermal fluctuations of the magnetization 
within individual interstellar grains when most interstellar 
Fe exists in a moderately ferromagnetic material (Draine 
& Lazarian 1998b). This last mechanism might be nearly 
independent of frequency in the 20-100 GHz range (see 
Fig. 7 of Draine & Lazarian (1998b) for Fea04 grains) 
avoiding the excess at 31.5 GHz described above. At least, 
one form of contribution — free-free or rotational emission 
or any other — is necessary to explain the correlations be- 
tween the anisotropies and far-infrared maps (Kogut et al. 
1996a; de Oliveira-Costa et al. 1997; de Oliveira-Costa et 
al. 1998). The role played by rotational emission with a 
spectral index a less than zero can be substituted by the 
free emission or others. It was shown that Galactic cirrus 
emission is high enough to explain observed anisotropies 
at 90 GHz (Fig. |J). Whether this level of anisotropy may 
be maintained for lower frequencies down to 50 GHz would 
be merely a question of fitting some parameter to other 
emissions. We can even fit the three or more emission 
types (free-free, dust rotational, dust thermal emission, 
magnetic dipole emission from dust grains,...) at the same 
time in order to get a nearly flat continuum spectrum be- 
tween 50 and 90 GHz, as well as to get a conspicuous 
increase for lower frequencies. 



where T is the antenna temperature, of the radiation re- 
ceived in a beam of FWHM=5 deg □. 

The aspect of the anisotropies is similar at all frequen- 
cies, and the widths of the peaks are similarO. However, it 
is usually claimed that anisotropies between 20 and 100 
GHz are predominantly cosmological while the other fre- 
quencies are dominated by the Galactic contribution. 

We also see this likeness in Fig. [| The correlation 
becomes null — first zero — around 30 deg at both Galac- 
tic and observed COBE-DMK. The first zero is related 
to the average angular size of the clouds when they are 
responsible for the anisotropies (Lopez-Corredoira et al. 
1998), and this coincidence implies that the average an- 
gular size of the supposed cosmological structures is ap- 
proximately equal to that of the nearby clouds giving rise 
to anisotropies. 

The linear size of giant clouds is normally between 20 
and 100 pc (Scoville & Sanders 1986) and can even reach 
sizes between 60 and 300 pc (Magnani, Blitz & Mundi 
1985; Heiles, Reach & Boo 1988). The projection of nearby 
structures such as these gives rise to inhomogeneities and 
irregular arcs extending between 10° and 50° on the sky 
(Low & Cutri 1994) in the form of infrared cirrus. In the 
range from 20 and 100 GHz, the Galactic contribution 
must also provide anisotropies with a first zero around 
20 or 30 deg because this is the mean angular size of the 
clouds, independent of the frequency (the first zero for syn- 
chrotron emission is somewhat less than 20 degrees as seen 
in Fig. ||a). This is precisely what is observed for the total 
anisotropies. Most authors claim that these anisotropies 
are cosmological rather than Galactic but the coincidence 

4 To get the equivalent antenna temperature obtained from 
the Tenerife Experiment with a FWHM=5° beam we do a 
convolution with a Gaussian response with a = 2.1°. 

5 There is a blank strip in IRAS 100 fim map, which inter- 
cepts S = 35° at a — 173°. Thence, the downward peak of 
Figure @ e) near this right ascension should not be considered. 
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Fig. 8. Second differences of antenna temperature in different 
surveys at S — 35°: a) 1.42 GHz from the Max Planck Institute 
for Radio Astronomy (Reich 1982; Reich & Reich 1986); b) 4.85 
GHz from NRAO (Condon, Broderick & Seielstad 1991; Con- 
don, Griffith & Wright 1993; Condon et al. 1994); c) 15 GHz 
from the Tenerife Experiment (Gutierrez et al. 1997), solid line, 
and 53 and 90 GHz from COBE-DMR (Bunn, Hoffman & Silk 
1996), (Figure taken from Gutierrez et al. 1997); d) 240 ^m 
(1250 GHz) from COBE-DIRBE (Boggess et al. 1992); e) 100 
^m (3000 GHz) from IRAS (Wheelock et al. 1994). 
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referred to here might be due to more than just chance. 
Though this does not prove anything, it is nevertheless a 
sign in favour of the Galactic predominance in MBRAs. 

3.3. Correlation among structures and the power spectrum 
at different frequencies 

Although anisotropics over the whole range of the electro- 
magnetic spectrum may be due to inhomogeneities in the 
Galactic distribution of dust and gas, this does not mean 
that the flux maxima and minima at different frequen- 
cies must occur at exactly the same coordinates. Several 
effects synchrotron, free-free, dust emission, etc. are re- 
sponsible for continuum emission but some effects predom- 
inate over others at different frequencies. Different effects 
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arise in different regions: synchrotron is higher where mag- 
netic fields are stronger, free-free radiation where a warm 
ionized medium is present and dust where the coldest tem- 
peratures are reached (Bennett et al. f 992). Peaks of dust 
emission are also to be observed at different positions with 
different frequencies since cold or warm dust are in dif- 
ferent locations; small particles — which are dominant at 
microwave frequencies (Draine & Lazarian 1998a) — and 
large particles — which are dominant in the far infrared 
(Greenberg & Li 1996) — may be distributed differently 
in the clouds, etc. Thus, we cannot expect uniformity in 
Galactic structures at different frequencies, i.e. an exact 
correlation for different frequencies. Such a nonuniformity 
is observed, for instance, by Davies, Watson & Gutierrez 
(1996, their Fig. 10). 

In any case, the correlation among different frequen- 
cies is not totally null. Comparison of the different plots 
in Fig. H shows a certain correlation. In Fig. 0, some 
cross-correlation is also observed at scales between 8 and 
30 deg, while there is an anticorrelation for less than 
8 deg. The microwave continuum in the range between 
14 and 90 GHz was also found to be correlated with 
100/im thermal emission from interstellar dust (Kogut 
et al. 1996a; de Oliveira-Costa et al. 1997; de Oliveira- 
Costa et al. 1998). This was interpreted as a good correla- 
tion existing between dust and free-free emission (Kogut 
et al. 1996a, 1996b); however, the correlations between 
Ha — normally a good tracer of free-free emission — with 
CMB and DIRBE maps are weak (Leitch et al. 1997; 
Kogut 1997). Leitch et al. (1997) alternatively proposed an 
anomalous bremsstrahlung emission from hot gas, but this 
was again inconsistent with the observed power radiated 
(Draine & Lazarian 1998a), so other kinds of emissions 
correlated with dust must be present. 



Hence, it must be concluded that correlation among 
different structures for different frequencies cannot be an 
argument either for proving or disproving the Galactic 
origin of the anisotropies, although some non-negligible 
Galactic contamination is necessary in order to explain 
these last correlations. 

The statistical distribution of the fluctuations given 
by the two-point correlation function, or its Fourier 
transform — the power spectrum — is also expected to vary 
at different frequencies. The coldest parts of the clouds 
produce dominant emission at lower frequencies so the 
cloud shapes vary when the frequency varies. Moreover, 
the dust power spectrum is very sensitive to the region of 
sky selected (Guarini, Melchiorri & Melchiorri 1995), the 
background subtraction and other details. As a matter of 
fact, different spectral indices have been obtained by dif- 
ferent authors for the dust emission: —3.5 < n < —2.5 
(Gautier et al. 1992), -2.5 < n < -1.5 (Melchiorri et 
al. 1996), n = -2.5 for |6| > 45° (Schlegel, Finkbeiner & 
Davis 1998), etc. 

At microwave frequencies, we cannot expect the same 
power spectra as in the far infrared. Some unknown 
changes are expected. According to this, power-spectrum 
or two-point correlation function shapes should not be 
used to determine whether the anisotropies come from 
the Galaxy or are cosmological. The fact that the angu- 
lar power spectra observed at high galactic latitudes by 
COBE-BIRBE are steeper than the COBE-DMR power 
spectrum should not be used, as in Wright (1998), to eval- 
uate the degree of Galactic contamination. The shape of 
the correlation function is also different from the obser- 
vational data in the microwave region in Fig. ^ there is 
a deficit of correlation at 9 ~ 10 deg and an excess of 
correlation at 9 ~ 25 deg. This may be, as has been said, 
because the shape of dust in far infrared was extrapolated 
in spite of the fact that variations were expected (§3.3), 
and that the cloud intensity fall-off from its centre is dif- 
ferent at 240 /im from that at microwave frequencies. 

3.3.1. The magnitude of the effects of the two-point cor- 
relation function shape variation with frequency 

The extrapolation of dust emission as a combination of 
three factors was carried out in (|J) under the assumption 
of a non-angular dependence of the first and second factor. 
This dependence would introduce much more complex cal- 
culations for which we do not have accurate enough data 
(it was not got sufficient information about the correla- 
tions at wavelengths longer than 240 yum) , although it can 
be proven that the order of magnitude of the previous 
calculations would not change. In any case, the approach 
taken here is better than any other prior to this paper and 
the results are more trustworthy. 

According to the hypothesis that clouds colder than 
the diffuse interstellar medium produce the anisotropies, 
the variation of the temperature within these clouds would 
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produce variations in the < TT > (9) shape. In a simple 
model of emission, in which the cloud flux is proportional 
to v 2 B(v, Tt), where Tt is the effective grain temperature, 
the antenna temperature T is proportional to B(v, T t ). T t 
is large enough to take the Rayleigh-Jeans approximation 
in the microwave range, so T cx T t . This implies that the 
variation of the antenna temperature in the microwave 
region is proportional to the variation of the effective grain 
temperatures. For the calculation of < TT > (6>), our 
approach contains a relative error 



S{< TT > (9)) _ 2 _AT t (0) 



(14) 



<TT> (6) " < T t > ' 

where AT f (6) is the mean variation of T t with respect its 
average value along two lines of sight separated by an an- 
gle 9. This is between 0.7 and 1.3 for a maximum relative 
variation of the effective grain temperature depending of 
the angle of 15%, which is quite a reasonable value (for 
instance, clouds within a range of temperatures between 
10 and 16 K with a mean temperature of 13 K produce 
a T t {6) between 13 K and 16 K, i.e. 14.5± 1.5, a 10%). 
This may justify the difference of the shape in the plots of 
Figure ^. In any case, the mean amplitude over the whole 
range is the one calculated above; there will be some an- 
gles in which the amplitude would be higher and others 
in which it would be lower although the Galactic contam- 
ination is within the order of magnitude of the observed 
anisotropics. 
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dence between the ratio of flux excesses and the filtering 
scale 9 (Fig. g): 



Ratio = (0.036 ± 0.024) + (0.0021 ± 0.0011) 9 



(15) 



Since this ratio as a function of the filtering scale is 
equal to 4 J f(9) dA, where A = ir9 2 and f(9) is the ratio 
as a function of the correlation angle, 



f(0) = (0.036 ± 0.024) + (0.0031 ± 0.0016) 9. 
On the other hand, 

T {ylB{ Vu T t {9))^ 



(16) 



f(0) = 



T 



100 /im 



T 



240 ^ r(viB{v 2 ,T t (9))^ 



2hi4 



e kT t(8) 



21u, 



(17) 



e fenW 



where v\ and are respectively the frequencies corre- 
sponding to 100 /im and 240 /im. From @ and @, we 
derive a linear dependence (correlation coefficient with a 
linear fit: 0.9958): 



T t {9) = 14.5 + 0.129 9. 



(18) 



This means that the range of temperatures for angles 
9 less than 30 degrees is within T t — 16.4 ± 1.9 K, i.e. a 
relative error of 11.5% and the mean temperature of the 
clouds would be (T t ) = 14.5 K (more or less in agreement 
with the calculations by Lagache et al. 1998). If we use 
equation (^), this leads to a maximum error of < TT > 
(9) of 23%. This may justify the difference of the shape 
in the plots of Figure ^, and it does not change the order 
of magnitude of our calculations. An exact calculation is 
not carried out since the uncertanties derived from Fig. 
U are too much and the use of a simple model of dust 
(v 2 B(v,T)) may be not totally correct to extrapolate T* 
from the given frequencies to the microwaves (for instance, 
very small particles may contribute less than a 10% of the 
total flux (Greenberg & Li 1996) which is negligible for our 
required accuracy). In any case, this was just to estimate 
the relative variation with 9 and not to calculate the exact 
value of T t . 



Fig. 9. Ratio of flux excesses (T*), i.e. after mean background 
subtraction, between the maps of IRAS 100/im and DIRBE 
240/i at |6| > 20° with different filtering scale 9. The solid line 
is 0.036 + 0.0021 9. 

From the IRAS 100/im and DIRBE 240/xm maps, we 
derive roughly (because the noise is quite high) a depen- 



3.4- MBRA galactic latitude dependence 

From Fig. pi it is observed that there is some positional 
dependence of the fluctuations. In this figure, the flux 
derivative is shown instead of the flux because it allows 
the fluctuations to be seen more clearly. The variation of 
the fluctuations at the Galactic poles is something differ- 
ent from that at intermediate latitudes. 

The density distribution of Galactic clouds is irregu- 
lar with a higher concentration in the plane and a fall-off 
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Fig. 10. Flux derivative versus Galactic latitude (6). a) 1.42 
GHz; b) 240 fim. 

240 um 



nates, and the column density follows another dependence 
with respect to b. The number of clouds is too small to pro- 
vide good statistics, so not a lot more can be said about 



this dependence. Figure 10 may show very slight trends 



but these are not too clear. 

If we do all the calculations of this paper in different 
regions, for instance: \b\ > 30°, |6| > 40° and |6| > 60°, 
the correlations expected would be somewhat lower. The 
anisotropies from dust are shown in Fig. [ll[ Instead of (Q) 
as first factor, the approximate fits derived from Reach 
et al. (1995) are (values for \b\ > 40° were calculated by 
interpolation): 



T dust (v,\b\ >30°)= ( 
xB(u,T= 17.70 K) + 
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which can be roughly represented by an exponential with 
a scale height around 75 pc (Scoville & Sanders 1986). 
Clouds are distributed very close indeed to the plane. As 
a matter of fact, nearly all clouds are less than 20 deg 
from the plane. The remaining clouds that we observe are 
a few local ones belonging to the plane and are very close 
to the Sun (Blitz 1991); they are distributed randomly. 
Any distribution in which the density depends only on z, 
the distance to the plane, should give a column density 
proportional to cosec{b), but in this case the irregularity 
makes the density depend also on the other two coordi- 
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T dust (v,\b\ > 60°) = (- 
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(21) 



And the values for the second factors due to dust ther- 
mal emission, calculated as above, are 3.5, 4.0 and 4.8, 
respectively, for |6| > 30°, |b| > 40° and \b\ > 60°. At 
90 GHz, rotational dust emission and synchrotron fluc- 
tuations are negligible, and those from the dust thermal 
emission, derived according to (|4|), predict a lower ampli- 
tude of the fluctuations at higher latitudes, as is shown in 
Fig. |j| 

The fluctuations, proportional to the square root of 
(TT), are about twice as high for \b\ > 20° than for 
|6| > 60° at a scale of around 8 = 5°. Therefore, a de- 
creasing fluctuation at higher latitudes is expected if the 
Galactic emission was the sole or at least an important 
contributor to the anisotropies. Perhaps the observed fac- 
tor between |6| > 60° and \b\ > 20° is not exactly two, 
since this quantity is subject to the errors of the calcu- 
lations of the first and second factor, but the order of 
magnitude is something of this order and there is likely to 
exist a decreasing correlation with Galactic latitude. 
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working. Until now, it cannot be said what exactly the b 
dependence in the anisotropics at 90 GHz is. 

Another con- 
spicuous dependence is that the anisotropics in the south- 
ern Galactic hemisphere are higher than in the northern 
Galactic hemisphere. This is still not a proof against the 
Galactic predominance of anisotropies at microwave fre- 
quencies because the COBE-DMR survey also observes 
higher anisotropies in the southern Galactic hemisphere 
(see Fig. 1 of Bennett et al. 1996)0. 

From this, the conclusion is that there is not a qualita- 
tive difference between the cosmological anisotropies and 
the cloud anisotropics. The position dependence of the 
MBRA is qualitatively similar in both microwaves and far 
infrared, so this may not be a proof neither for the cos- 
mological origin nor for the Galactic origin. 



Fig. 11. Angular correlation function at 240 /im. 
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Fig. 12. Predicted Galactic two-point correlations at 90 GHz 
at different zones. 

In fact, MBRAs at microwave frequencies show a 
slightly decreasing dependence on \b\ with the fluctua- 
tions (Smoot et al. 1992, Fig. 2). They show that the 
correlations for \b\ > 30° and \b\ > 40° are nearly a half 
that for \b\ > 20° for angles 9 less than 10°, although 
their first-year COBE-DMR data are quite noisyB. Data 
showing this dependence more accurately are still awaited, 
specially when next experiments (PLANCK or MAP) be 

6 I do not undesrtand why these correlations were not calcu- 
lated again with the four-year COBE-DMR data but only for 
161 > 20°. 



3.5. Conclusions about possible Galactic predominance in 
MBRAs anisotropies 

The conclusion is that under some particular, but not im- 
possible, conditions, all the microwave background radi- 
ation anisotropies may be due to Galactic foregrounds 
rather than cosmological in origin. There are no argu- 
ments yet to exclude this possibility although this is not 
yet proved and the question remains open. A testable pre- 
diction of such a case would be that the amplitude of the 
fluctuations for |6| > 60° would be about the half of that 
for \b\ > 20° for angles around 5 degrees. 

The implications of such a question are extremely im- 
portant, not merely for refining some quantity or other or 
for making certain corrections to get an accurate result 
for an individual parameter, but because it would result 
in a different qualitative description of the Universe. The 
implications for inflation theories or the formation of the 
large-scale structure of the Universe would be enormous, 
and our ideas regarding such formations would change 
completely. Hence, I think studying the influence of the 
Galaxy is a valuable exercise, in order to avoid the haz- 
arding of cosmological theories based on cumulative errors, 
in which this paper claim to be still an open question. 

4. Conclusions 

The following main conclusions may be drawn from this 
paper: 

— The extrapolation of anisotropies following the mean 
dust emission is a bad approximation since it does not 
take into account the growing contrast of colder clouds 
in the background of the diffuse interstellar medium. 

7 I do not understand why this fact is not commented on in 
that paper, since the figure shows the difference of both hemi- 
spheres quite clearly. This difference is attributed to Galactic 
contamination, but may not Galactic contamination be respon- 
sible 100% of the anisotropies? 



14 



Galactic microwave anisotropics 



By considering this effect, it is found that dust ther- 
mal emission anisotropics are higher than expected by 
other authors, and that their amplitude is comparable 
to the observational data at 90 GHz. 

— Our ignorance of the different emission mechanisms 
around 50 GHz (free-free, dust rotational emission, 
magnetic dipole emission from dust grains) do not al- 
low the firm conclusion that anisotropics due to Galac- 
tic emission are not frequency dependent but this pos- 
sibility remains open. 

— If Galaxy-induced anisotropics are not responsible for 
the totality of the correlations they would at least be a 
non-negligible part of them, so untrustworthy cosmo- 
logical conclusions could be reached from microwave 
background radiation anisotropy analysis unless pos- 
sible Galactic emission processes are correctly sub- 
tracted. 

— If the Galaxy-induced anisotropics made up the total 
correlations at microwave frequencies, then inflation, 
models of Galaxy formation and many parts of the 
standard cosmology would be wrong. This is not im- 
possible, though there is no firm evidence either for or 
against it as yet. 
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